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CHAPTER 3 
a 
Wave Loading 
Introduction 
1 The loads caused by waves on surface ships are by far the most important in normal 


service and will dominate the design of the primary structure. The presence of a hull in the sea 
will result in a general lack of correspondence between mass distribution (m(x)) and buoyancy 
distribution (b(x)) along the length. If a ship was able to maintain an equilibrium position on a 
fixed wave profile then the total mass and mass moment must be equivalent to the buoyancy and 
buoyancy moment, that is: 


(3.1) a 
J (m(x) — b(x)) dx = 0 


(3.2) ae 
'- x (m(x) — b(x)) dx = 0 


where L is the waterline length and x is measured forward from amidships. Then at any given 
section x; there will be a bending moment given by: 


(3.3) x, 
f (c+ L/2) (m(x) — b(x)) dx = M, 
-L/2 


ns for a given sea surface profile is known as quasi-static balance, 
he motion of the ship this equilibrium is, of course, never achieved 
water. Consequently, the effects of translational and rotational 
accelerations, as well as the instantaneous lack of static balance, need to be added to the quasi- 
static calculation. This is a complex task and is discussed in detail later in this chapter, although 
the general problem has been known to naval architects for many years. 


The solution of these equatio 
Due to the dynamic effects of t 
except with the hull at rest in still 


a special case insofar as there are no dynamic 
tic calculation. Although, therefore, not strictly 
under the discussion of wave action to save 


- The bending moment in still water is 
considerations and it can be found by determinis 
due to wave action, it will nevertheless be covered 
later repetition of what is a standard procedure. 


ave loading has to be dealt with in a linear manner because 


3... The fu ic analysis of W 
¢ full dynamic analy table with currently available computational tools. 


the complete non-linear equations are intrac 
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even on iron hulls. Rankine also gives tab 


The assumption of linearity means that the hull is taken to be wall-sided to the ful] amplitude o 


the wave: this is, at any section the buoyancy per unit length 1s directly proportional to a 
elevation. Additionally, the wave surface is assumed to be the summation of a large number fe 
sine waves. with the calculation being carried out for single waves comprising a wave Spectrum 
and the individual responses being summed over all the sine waves within the spectrum. The 
effect of these assumptions is that the derived bending moment varies linearly with wave height 
for a given speed and heading, and there is no difference between sagging and hogging responses 
Loads are also overestimated because the hull is assumed to be ‘infinitely wall-sideg: and 
buoyancy is added or removed even when the wave profile is above the deck or beneath the keel. 
Ship motions in relatively benign sea states can be estimated with reasonable accuracy, py 
structural responses in the high sea states of interest to the structural designer are of doubtfy) 
validity. (Note that wave height is defined as distance from trough to crest, while amplitude is the 
maximum elevation from the still water line.) 


4 From the discussion above and in Chapter 2 it is clear that the question of the true 
magnitude of wave loading on a hull in a given lifetime is a very complex one. It has been 
approached with increasing sophistication over the years as greater numerical computing power 
has become available, but to provide a thorough understanding of the current position it is 
worthwhile to present it in a historical perspective. 


History of Wave Bending Prediction 


5 The first classical publication on Naval Architecture in the modern sense is by Rankine 
(1866) entitled “Shipbuilding, Theoretical and Practical’ and which foreshadows much of the 
developments of the next century. He is the first to define rigorously such terms as stress and 
strain, and also includes clear descriptions of the strength and mechanics of beams and trusses. 
His exposition of bending moment and shear force calculation is surprisingly modern and by 
dividing the ship into two halves for the estimation of bending moment amidships, he provides 
the basis for the much later ‘Murray’s method’ (paragraph 15), 


6 Rankine also presents the bending moment prediction formula: 
(3.4) M = ALIC 


later attributed to John (1874), where A is dis 
20 in hogging as an upper bound on M for 
mass distribution. Sagging at that time is alw 
ships of the day, and deck buckling was app 


placement and C is a constant suggested as equal to 
rectangular section hulls with an assumed standard 
ays dismissed as being less than hogging in the short 
arently not a problem as most decks were still of wood 


ular methods for obtaining the position of the neutral 


axis and the second moment of area. th idi 
| , thus providing all] esses ID 
the hull as a simple beam. P g all that is necessary for calculating str 


7 Reed (1871) goes further and 


res | 3 ss tas 
force and bending moment curves ie ents a number of examples of graphically dert ved sh 


various warships, both in still water and balanced on waves 
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In fact Reed appears to be first to us : 
St to use balancin a wave crest or in a tre st 
£ On a wave crest or in a trough as explicit worst 


cases for wave loading. He also presents an intuiti 
aa s | sents an intuitive argument for waves equ: 
ship giving the worst load. § es equal to the length of the 


8 The question of waves and the effect of dynamic loading was being pursued by a number 
oR okie toads ae late 19th Century, notably Smith (1883) and Read (1890) who present 
methods for allowing for the effects on water pressure (and therefore on bending moment) of 
particle circulation in waves, and of pitch and heave. These methods were developed with 
sophisticated graphical solutions by Kriloff (1898) and Alexander (1911). However, Robb (1918) 
shows that the extreme hogging curve for quasi-static balance calculation covers all other 
conditions and that the Smith and Read corrections have insignificant effect compared with the 
effect of changes in the mass distribution. Thereafter the use of corrections, especially that due to 
Read, becomes progressively rarer, although the Smith correction is once again coming into 
prominence as it may be significant in the extreme waves now used in design load prediction as 
discussed by Clarke (1986). 


9 Continuing with the development of the design process, John (1874) took the method as 
described by Rankine (1866) and, by applying it to a number of ships, deduced the formula: 


(3.5) M = AL/35 


which subsequently became known as the ‘Standard Bending Moment Formula’ or more lately 
‘John’s Formula’. The particular aim of this formula was to prevent the overstressing of large 
ships by relating bending moment to length as well as displacement. The Lloyd's Rules which 
were current at that time required scantlings to be directly proportional to displacement with the 
result that the factor of safety against yield fell from about 5 for a 100 ton ship to 2 for a 4000 ton 
vessel. It should be noted that John intended the formula specifically for merchant ships. There 
was no agreed constant for warships until the work of Champness (1922) confirmed the value of 


35 for hogging loads but gave 29 for sagging. 


gave an estimate of bending moment, there 10While the 


standard formula gave an estimate of bending moment, there was still a need for an agreed 
acceptable stress in the deck and keel, and at later stages of the design where static balance was 
used to estimate the bending moment, there was also aneed for an agreed wave height. On the first 
point there was no common consensus at the time and there appear to be no recorded values of 
acceptable stresses in any design documentation since. However, the notes provided for students 
of Naval Architecture at the Royal Naval College, Greenwich in the early 1960s require that the 
maximum deck compressive stress in a mild steel frigate (of about I 00 m length) should be not 
more than 5.5 tons/sq.in (85 MPa) and at the sketch design (first iteration) stage 4.5 tons/sq.in (70 
MPa) should be the limit. The notes also state that higher stresses may be acceptable in higher 
strength steels, but that in compression the advantage may be limited by buckling effects. 
Although the values are quoted for a frigate, clearly material properties are independent of ship 
size (except insofar as properties vary somewhat with material thickness - see Chapter 5) and 
these values are a good guide for any designer of a large steel structure at the beginning of the 


- design process. hee 


10 While the standard formula 
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13 The initial design process was then to estimate 
formula, equation 3.5, and using an allowable stress 
(Z) at amidships from the beam formula: 


the bending moment from the standard 
(G) to derive the necessary section modulus 


(3.6) Z=MIo 


A section could then be synthesised by trial and error to match the required modulus, and then 
scaled lengthwise using experience to give the scantlings for the rest of the hull. The structural 
weight could then be deduced and added to all the other ship weights which had been estimated by 
comparable means, giving a complete ship mass distribution. A static balance was then carried 
out to provide a more accurate estimate of bending moment and the process repeated progress- 
ively to refine the scantlings. 


14 A variation on this method was proposed by Biles (1908) in which he suggested that for 
merchant ships a standard mass distribution could be used. That was a trapezium with the same 
mass and centre of gravity as that required for the ship and which became known as ‘Biles’ 
Coffin’ (Figure 3.2). It was used extensively to avoid use of the standard formula and to enable 
entry to the static balance calculation at the start. For warships with a much more uneven mass 
distribution Biles’ method was not appropriate and standard formulae continued to be used. A 
method for estimating mass distribution for warships based on a constant volumetric density has 
been proposed recently by Chalmers (1982 and 1988(b)) and will be presented in more detail in 
Chapter 6. 





L/3 L/3 L/3 


Figure 3.2 Biles’ Coffin Mass Distribution 


15 Historically, perhaps one of the most important milestones came with Murray CI 946-7) 
who collected a large volume of data and suggested a much more refined method for estimating 
bending moment in which moments due to mass, still water buoyancy and wave buoyancy are 
each calculated or estimated independently. The wave bending moment, defined as the difference 
between the still water condition and that with the ship on a wave, is proportional to L3B (Bis the 
beam amidships), the still water bending moment is proportional to L (as John s formula) but is 
also linearly related to block coefficient (C,,), and the moment due to mass is calculated from an 
assumed known mass distribution. An estimated distribution such as Biles’ Coffin may be used if 


the real distribution is not known, but Murray does not recommend it. 


he most useful for merchant ship design but unfortunately, 
lex hull form and mass distribution, the method has never 
s. Baker in his discussion to Murray’s paper (1946-7) 


16 Murray’s method is still one of t 
presumably due to the much more comp! 
been successfully developed for warship 
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tries to fit a number of warship designs in Murray’s theory but fails to get any reasonable 
correlation. He suggests that relating the still water bending moment to prismatic coefficient (¢ 
rather than block coefficient might be more reliable for warships, but the idea has Never been 
followed up. Murray’s curves for still water bending moment do not extend toa low enough Value 
of C,, for warships and simple extrapolation does not seem to give sensible answers: Lioyd’s 
rules, where they quote the method, say that it should not be used for values of C,, less than 0.6, 


17 At about that time the use of statistics began to be proposed for the prediction of design 
bending moments. A number of publications appeared through the 1950s and 1960s of Which 
Turnbull (1953), Lewis (1957), Telfer (1960), Yuille (1963), Ward (1963), Murray (1 963), Moor 
(1967) and Murdey (1972) are all worthy of study as they increase understanding of the subject. 
although none have had any lasting effect on commonly used design processes. The text books 
written in the same period, in particular Barabanov (1967), Muckle (1967), Harvey Evans (1975) 
and Comstock (1977) all still advocate the use of John’s (1874) method or Murray’s (1946-7 } 
method for initial estimation of bending moment; the use of the simple beam formula and an 


allowable stress to give an acceptable section modulus and Static balance on, usually, a 1.1VL 
(0.6VL metric) wave for subsequent iterations. 


18 In recent years a considerable amount of work has been put into the development of 
numerical methods for derivation of design bending moments so as to be able to allow for actual 
wave form, relative ship motions in waves and the dynamic response of the ship, rather than a 
notional quasi-static condition. These methods are known generically as strip theory; the ship is 
divided into vertical strips and the linear response of each strip to a sinusoidal wave excitation is 
calculated while satisfying the compatibility conditions between each Strip. For a complex or 
irregular sea linear superposition is used, as mentioned in 
achieve the desired sea state. Various methods are availa 
due to Gerritsma et al (1967), Salvesen et al (1970) and Vugts (1969), and a number of equivalent 
computer programs have been written to solve the numerical equations, for example SCORES 


the ship to be rigid beam: more recently 


ble, the most commonly used ones being 


linear analysis and extreme effects are not well m il 
strip theory is discussed later (paragraphs 3] 
_ Into the non-linear regime, for example J 

_ academic interest, | 


-49). Efforts have been made to extend strip theory 
ensen and Pedersen (1979), but they remain mainly of 


__ Current Available Prediction Methods 
eae ns wlio the foregoing historical survey it is clear that the modern designer has to choose 
Ser aceng ecg etnctbeen meet en pass 

BRU pases nee itly OF the dynamic effects of ship motion in waves, or embracing more 
sae rte soi esi based analyses. which allow dynamic response but lose the effects of non 
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ore close somp; ; 
osely and to compare the results of analyses with 


vice. These latter data are, after all, the only hard 
Out what it actually happening at sea, . 


measurements taken from ships in ser 
information available to the designer ab 


ct le etl Static balance calculations in relation to the need or otherwise to allow for 
’ investigated by Robb (1918) who showed that the standard calculati 
without any allowance for heave and pitch motions was conservative, and that snitch sainets 
reduced the estimated bending moments. However, his work related fosthie relatively low, hnagh 
severe, waves around the L/20 height at which hogging bending moments are very citer 46 
sagging ones; indeed, for the ships of that time hogging was seen as the worst case. In the extreme 
conditions which relate to the Ultimate Limit State, that is the point at which the hull will fail 
irreversibly (contrasting with the Serviceability Limit State when the capability of the hull starts 
to be degraded, but may still be safe), the worst condition is Sagging as illustrated in Figure 3.3. 
This shows, for static balance on increasing wave heights, that the hogging bending moment falls 
below the linear relationship which would be expected for a wall-sided hull, at a much lower 
wave height than the sagging response. Measurements from sea presented by Brown (1977(b) and 
1981) and Clarke (1985) and model tests presented by Brown (1978) and Lloyd et al (1980) 
confirm this conclusion for real conditions. It is clearly the extreme conditions which will damage 
or sink ships, and for the designer to have an understanding of the margins he has between design 
loads and failure so that he can obtain the optimum performance for his structure, it is necessary to 
investigate whether Robb’s conclusions are equally valid in the extreme sagging state. 


21 A study of Robb’s (1918) paper shows that the sagging case, although of lesser magnitude 
than hogging, is significantly affected by heave and pitch. Indeed, Robb himself notes the fact and 
remarks that for ships of high speed form such as liners or ferries (and therefore including 
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~ Figure 3.3. Comparison Between Extreme Hogging and Sagging Bending Moments 
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warships) his conclusions may not be valid, especially if the vessel has heavy machiner, 
amidships increasing the still water load. Unfortunately, as noted earlier (paragraphs 3 and jg) 
strip theory gives no help in such comparisons as the hull form is assumed wall-sided and hoggin 9 
and sagging responses are identical. In the absence of any non-linear theory, which would in am, 
case be numerically extremely complex and need to operate in the time domain (as opposed to the 
frequency domain for linear strip theory) the best solution to the problem is to compare static 
balance results with measurements taken from ships at sea, and to study the correlation. 


TABLE 3.1 


Comparison Between Expected Maximum Bending Moment Based on Measured Strains and Static Balance on a 
Wave of Height H = 8 m 





BM from Static Balance Ratio 
(MNm) Measured/Calc 

Ship Class H/L Hog Sag Hog Sag 
Type 12 (Rothesay) 0.073 122 180 1.02 0.97 
Type 81 (Tribal) 0.075 111 160 0.98 1.14 
Leander (N Beam) 0.073 123 177 1.12 1.17 
Leander (B Beam) 0.073 127 194 0.97 1.1] 
County 0.052 325 425 1.00 1.05 
Type 21 0.073 139 221 1.08 1.1] 
Type 42 (Batch 1) 0.067 165 264 0.81 0.90 
Type 22 (Batch 1) 0.064 188 282 0.90 0.97 
CVSA 0.042 2200* 0.96* 
eee eee eee es, wy RUSS lhe eee ee MI 
Average (not including CVSA) 0.98 1.05 
* Total range +0.10 +0.10 
— 
22 


The only published work in this area is due to Clarke (1985 and 1986) who has taken strain 
measurements from a number of ships in the length range 100 m to 200 m and derived from them 
the expected bending moments (the values that are expected to be exceeded at least once) in a ship 
life of 3 x 10’ wave encounters (see Chapter 2 for a discussion of ship life). These results have 
been compared with static balance calculations on an 8 m wave without the Smith correction (for 
numerical convenience but see paragraph 8) and are listed in Table 3.1 which demonstrates good 
correlation for both hogging and sagging. The derivation of the measured bending moments is 
explained in the references, but briefly it is based on Strains measured on a number of shi ps of the 
classes listed and plotted in a Gumbel (1958) form. The design value which is compared with the 
static balance calculation is first converted to bending moment from strain by aed beam theory 
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this strip theory calculation is achieved by scaling from static balance results but even so the 
sagging prediction does not correlate as well as the results presented in Table 3./ 


23 However, strip theory does take account of phenomena which cannot be explicitly 
included in static balance bending moment predictions and it may therefore be useful for 
providing an insight into unusual happenings resulting from dynamic responses which are related 
to non-extreme and so relatively linear motions, for example the effects of different modes of 
operation on fatigue loading. while it may appear therefore that a version of static balance iS 
currently the most accurate prediction method, neither method is particularly good at predicting 
extreme load values, and strip theory undoubtedly has more scope for development. Both 
methods are consequently presented in the following paragraphs, in each case with a procedure 
for extrapolating the result to an extreme condition appropriate to the design value, that is the 
ultimate limit state. 


Design Bending Moment by Static Balance 


24 Static balance theory is straightforward and may be found in many Naval Architecture text 


books, for example Rawson and Tupper (1976) or Comstock (1967), and in any case the method 
is now embodied in well established and validated computer programs. It will be assumed, 
therefore, that the designer is able to carry out a static balance calculation on any height of 
trochoidal wave for a conventional monohull structure. The aim then is to achieve a resulting 
bending moment which reflects an expected probability of exceedance in a given number of wave 
encounters. As discussed above (paragraph 22) for warships in the length range 100 to 200 
metres, the expected maximum bending moment in 3 x 107 wave encounters for British naval 
areas of operations can be approximated to by balancing the ship on an 8 metre wave of the same 
length as the ship. However, this is an expected value, not a design value, and there are other 
uncertainties implied by the standard deviations in Table 3.1. Furthermore, the strain measure- 
ments from which the 8 m wave height was deduced were taken from ships with a worldwide 
operating history. The conditions in the North Atlantic are such that ships of around frigate and 
destroyer size meet waves of their own length with a greater frequency than in most other sea 
areas. There is evidence that ships in this length range (100 to 150 metres) could, for the same 
probability of exceedance, be subjected to lifetime design bending moments up to 15% higher 


than the worldwide average if operating solely in the North Atlantic. 


ocedure for estimating design bending moments is probably best illustrated 
ner can insert his own parameter values where those differ from the 
following. The only critical parameter which needs to be decided and on which most of the rest 
depends is the required probability of exceedance of the design bending moment in a ship life. A 
study of failure rates in other engineering disciplines, see for example Brown and Chalmers 


(1990), particularly those associated with large steel structures, suggests that a figure of the order 


of 1% is about average, and that is what will be taken here as an example. It is therefore a question 


of investigating the uncertainties associated with the 8 m static wave balance, and then 
extrapolating the result from the expected value to that which has a 1% probability of exceedance 
in 3 x 107 wave encounters. For other numbers of wave encounters see paragraph 30a. 


25 ~ Asuggested pr 
by an example and the desig 
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Figure 3.5 Probability Density Distribution of Maximum Hull Girder Bending Moment 


the 1% point. As is described by Clarke (1986), the actual measurements above the expected 
point, insofar as they exist, tend to drop below the straight line due to non-linearities in the 
extreme response. The effect is to change the distribution about the expected value from a skewed 
one to an approximation to a normal one as shown in Figure 3.5. The standard deviation of this 
‘quasi-normal’ distribution is 0.15. Combining this deviation with that above, assuming they are 
independent, gives a value of 0.18. Note that, as explained in paragraph 24, a worldwide average 
design bending moment is still being estimated. A further partial factor will be needed if a 
significant proportion of the ship’s life is to be spent in the North Atlantic. 


28 To achieve a 1% probability of exceedance 2.33 standard deviations are needed and so 
combining this with the bias gives a factor 1.12 + (2.33 x 0.18) = 1.54. This factor only applies to 
the wave bending component of the static balance results because the measurements on which the 
correlation is based are zeroed at the still water condition. It may be assumed that the still water 
bending moment can be calculated directly with no significant error, although there will clearly be 
some scatter due to uncertainties in the displacement and mass distribution and the fact that there 
will be some variation with ship speed, as shown by Brown (1978) and Lloyd et al (1984). 
however, still water bending moment is only about 10% of the wave bending moment and so any 
inaccuracies will have relatively little effect. Nevertheless, the designer should try to estimate the 
worst still water bending moment to which the ship is likely to be subjected over its lite. 


29 Design bending moments may then be deduced from the following formulae: 
(3.7) Mj, = Moy + 1-54 (M, — Mow) 

(3.8) My, = Moy + 1.54 Ma - M,,,) 

where M,., M,,, are the design sagging and hogging bending moments; 


M,,M,, _ are the sagging and hogging bending moments from static balance on an 8 m 
s? 


wave; | oh 
 M, is the still water bending moment making proper allowance for its sign 
Ww - ‘ 
(usually +ve for hogging and —ve for sagging). 
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These expressions can be used to derive bending moments at all points along the length of the 
ship. However, to allow for slamming effects the maximum value at amidships should be 
extended forward by 15% of the ship length and then reduced linearly to zero at the fore 


perpendicular as recommended by Clarke (1986). 


30 This process has been validated against a number of ships designed by more traditional 
methods where the bending moments used at the time of the design were found from balance on 
an L/20 wave and a simple stress factor used as an acceptance criterion. The ultimate strength 
values to which the design bending moment is compared, and which are intended to relate to the 
ultimate limit state (see paragraph 20), are defined and explained in Chapter 7. The overall 
implication is that ships designed to a load with a 1% probability of exceedance by the above 
method will be slightly stronger than previous successful practice. 


31 The number of 3 x 107 wave encounters was chosen to reflect a typical RN warship life as 
indicated in Chapter 2 (paragraph 10). If however the ship is to be designed for a longer or more 
intensive life then 3 x 107 wave encounters will be insufficient. The logically correct way to allow 
for more wave encounters would be to increase the static balance wave height, but to achieve this 
would mean repeating the previous work with different parameters. A simpler and sufficiently 
accurate method is to increase the 1.54 factor in equations 3.7 and 3.8. This can be done by using 
the process of Clarke (1985) to derive the 1% exceedance probability for longer lives and then 
using the slopes of the lines in Fig. 3.4 to provide a ratio for longer lives. Using this derived ratio 
the factor to be used in equations 3.7 and 3.8 is then as shown in Table 3.2. 


TABLE 3.2 
Bending Moment Factor for Higher Numbers of Wave Encounters 


No. of wave encounters Factor in equations 3.7 & 3.8 


3 x 107 1.54 
4x 10? Lo 
5 x 107 1.59 
7x 107 1.63 
10 x 107 1.67 





It is possible that these factors are a little optimistic due to the greater uncertainty implicit in the 
further extrapolation, but this should be offset by the known trend for stresses to fall below the 
linear extrapolation at lower probabilities. 


Strip Theory 


32 There is no particularly clear explanation of strip theory available and probably only that 
due to Bishop and Price (1979) is readily accessible. Consequently, a brief summary of the 
method will be given here. As the loading response of the hull is a function of both the inertial and 
fluid forces, and these are functions of the ship motion, it is first necessary to solve the rigid body 
equations of motion. . Foe | | 
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Figure 3.6 Strip Theory - Axis System and Nomenclature and Sign Convention 


33. Bor the axis system shown in Figure 3.6 the following are the equations of motion for 
heave and pitch, neglecting, for simplicity of presentation, the damping terms z and ©: 


(3.9) \; A = BA 
I,8|° LM(1) 
where A= mass displacement 


L= pitch inertia 


F(t), M(t) represent the wave induced force and moment acting on the ship and are 
determined by integration of the wave induced fluid loading oF (x,t) as 
indicated in Figure 3.6. Note that fluid damping effects are included on the load 
side of the equation, as stated in paragraph 35, to simplify the solution so that 


(3.10) F(t) ef LD dx and M(t) = | FD 5 ae 
hax Ox 


34 In elementary strip theory ‘t is assumed that the wave induced fluid loading for a strip at x 
is a function of the relative vertical displacement z,(x,t) such that 


(3.11) z,(x,t) = 2(%t) — C(~0) 
where z is the absolute displacement and G is the wave elevation at x (+ve up). 


35 _ It is further assumed that the fluid loading is made up of three terms, the fluid added mass, 
fluid damping and the effect of buoyancy. The first two terms are functions of strip shape, 
encounter frequency and relative vertical displacement and details of their computation can be 
found. in Bishop and:Price (1979). The third term is simply pgB(x)z,(x,t) where B(x) is the 
(constant) beam at X. 


36 The fluid force may then be re-expressed as functions of z and C so that 
3.12) OF _oJ OK 
: ax ox ox 
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It also follows that the 


area under the distribution curve is equ 
that 


al to the mean square value m,, such 


(3.24) 62 ee 
m, = | ®..(@) dw = i 72 aa 
0 


r=0 


43 From lengthy observations, for example those due to Hogben et al (1986), the sh 
wave spectra can be correlated to two parameters, the si AE, 

characteristic period (T). H,,3 is the mean of the highest third of all observed waves and is equal to 
3.8Vm, (a little lower than 4.0Vm, given by a true Rayleigh distribution; see International Ship 
Structures Congress (ISSC) Report 1.1 (1982) for a more detailed explanation). ®,,(@) is 
therefore a function of H,,3 and T, and can be found from a number of empirical formulae, the one 


most commonly used in structural design being the ISSC, ITTC (International Towing Tank 
Conference) or Bretschneider two parameter formula, thus 


ape of 
gnificant wave height (H,,,) and a 


3.29 A —B 
sera D_.(@) = as °XP ka 


where A = (173 H?,,)/T} and B = 691/T} 


There is also available the JONSWAP (Joint North Sea Wave Project) spectrum for limited fetch 
conditions, most appropriate for ships designed exclusively for restricted waters. This spectrum 1s 
defined as 


(3.26) ® (@) = 0.658 C P,(0) 


where C = 3.37, g = exp [(-1/y2) {(@T,,/27) - 1 }2] and y= 0.07 if a on - is 0.09 a ae 
The characteristic period T can be described ina variety of ways and for if te re . ae 
formulae for A and B the characteristic mean period T 1s assumed but in t e ee 
the modal period T,, is used. Wave statistics may be quoted against a ee is < te Hie 
which may be aunverted to T, or T, as follows taken from ISSC Report I. 


(1981): 


Modal period T, (or T,) T,= Oat a} 2 5 wi st 
Zero upcrossing period T, T, = ree A 5 ze 5 ee Ne 
Significant period Tj/3 T, = 9. 1/3 ° 


taken from seas in many 
de scatter of measurements 
i averaged over a WI! acho 
see ae ote if Sento is to be created fora known sea ne sine eng po 
sath Hep P esaieiits, However, of all the different periods quote at fou hue 
wae aneniey (sometimes called T,) and therefore that which will be mos 
0 the visual period (s 


i iod. 
data has been collected visually, 1s the zero upcrossing pero 
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(3,27) PD (OX) = Dp (a), = (X) 


4 function describing the spread of the waves about the dominant direction Ko 
i FOr 


where = (0) 8 | | 
etrum quoted in pal 
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, 


s AO 2 , 
(3.238) S(y)= : cosy 


is recommended, Further discussion on wave spreading can be found in Bishop and Price ( 1979), 
45 For any particular combination of speed, heading and encounter frequency the response of 
the ship can be derived from the wave input through the Response Amplitude Operator, which 
is the square of the Transfer Function (H,,(V,x.,). These response amplitude operators may be 
derived theoretically as described above for regular seas and for a linearised hull form, or may be 
found from model tests. A good description of the process is given by Comstock (1967). Thus, for 


a given regular wave, the response 
(3.29) R, (V,x.0,)2 =| A, (V.x,0,)17 a? 
and it follows that the response spectrum is given by 


(3.30) ®, (V.x,0,) =| H,. (V,x,,)? ®., (@,) 2 


The response spectrum may be in terms of motions (heave, pitch) or, for structural purposes in 
terms of shear force and bending moment, depending on which response amplitude operators are 
used. The total energy in the response spectrum, m,,, can be found by integration as follows and is 
then used in the statistics of the response. 


(3.31) oo Tt 
mor V)= FF ®,, (Vitae) dx de, 
Ome 


Tt 


hich 


46 If the bending moments experienced by a ship over a short period of tme during W ‘a 


conditions remain sensibly constant are measured at random or at equal time intervals, 
Statistical distribution of the measurements will approximate to a Gaussian (normal) distributio” 
cre for most structural purposes the statistical distribution of the maxima and minima Is ° 
soa epee = a narrow banded spectrum (which is usual for short term response a) 
this form of distri urning values follow a Rayleigh distribution (see Price and Bishop ( 986): 
This bias astnuton me confirmed by measurements from sea as discussed by Clarke - ich 
can be sh as Passa nis because the distribution is a function of a single paramere 
shown to be equal to the integral of the response energy spectrum, ine 
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response m,,. Thus the probability 
measured over a short period of time 


density function for the m 
in which conditions remain 


? _R2 
(3.32) P(B) = = exp =] 


-_ > — 
m <M, 


or 


aximum bending moment B. 
constant, 18 given by 


The probability distribution function. 


. or cumulative probability, of the maximum bending 
moment equalling or exceeding a given v 


alue of B in any one cycle under the specified conditions 


is given by 

(3.33) P (B> B)=exp (~B?/2m.,,) 

It can also be shown that the probability of exceeding the response B at least once in N cyeles is 
given by 

(3.34) Py (B > B) = 1 — exp {—N exp (—B2/2m.,,.)} 

47 


The short term statistics described above only remain valid as long as conditions remain 
steady. At sea this is unlikely to be for more than about an hour and so it is necessary to take 
account of the probability of meeting every combination of circumstances when using the process 
to estimate design loads. To begin with the probability of meeting all possible combinations of 
wave height and period (P(H,T)) in the two parameter spectrum will need to be established. It will 
then be necessary to assess the proportions of operational time spent at different speeds and 
headings with respect to the wave conditions, noting that courses and speeds will be constrained 
by severe weather and therefore their probabilities will be functions of wave height. Values of 
P(H,T) depend on the areas in which the ship will operate and data can be extracted from 
references such as Hogben et al (1986). The other two probabilities are less easily determined and 
perhaps the only source of data would be ships’ logs, leading to considerable uncertainty over the 
reliability of the process. 


48 = Assuming that such information is available, the probability that the peak bending 
moment B in any one cycle, which can now occur in any of a wide range of conditions, exceeding 
any given value B is given by 


(3.35) me 


P(B > B)=S SSS POA). PQ! H) . PCT) exp Ewa al dy dV dT dH 


This formula can be simplified somewhat if it is remembered that it has been assumed that 
responses are directly proportional to wave height. Nevertheless, the numerical computation of 
the probability is tedious and requires the input of a very large quantity of data, some of which 
may be of doubtful ancestry, particularly that relating to operational scenarios. The probability of 

- @ particular bending moment being exceeded in N cycles can be found from an extension of 
equation 3.33, | 


ie x elt aims 


=. 


_— 


49 si re of bending moment calculated in this way can be compared for particular 
| _ probabilities of exceedance with values measured at sea and plotted in the Gumbel form, Figure 
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3.4. However, the calculated loads are always significantly greater than those derived ¢ 
measurements, even though slamming is not included in the theory (see paragraph 5 
because of the assumption of geometric linearity, use of simple beam theory, 
assumptions and uncertainty as to the ship condition during strain measurements, an 
necessary to deduce some form of factor to use to achieve realistic design loads. 


TOM Strain 
1), Mainly 
sea State 
d so It js 


Allowance for Hull Flexibility 


50 _—_— All of the above theory assumes that the hull is acting as a rigid beam. The hull is of course 
in reality flexible, and it is possible to include terms in the equation of motion for each Strip (3.18) 
which take account of the shear and bending stiffness of the strip and relate shear force anq 
bending moment to the applied forces (Bishop and Price (1979)). Such analysis will yield the 
flexible modal responses as well as the rigid body responses, and if there are any distortional 
resonances in the frequency range of interest there will be further bending moments and shear 
forces to add to those due to direct wave motion. Under normal conditions the flexible responses 
are sensibly the same as the rigid body responses, but if slam or explosion induced whipping is 
included then the flexible analysis is essential for a complete understanding of the phenomena. 
Unfortunately, however, the responses of interest to the structural designer are well into the non- 
linear range and so current theory can only give a very approximate guide to likely loads or 
motions and should be used with great caution, but there are occasional incidents, for example 
those described by Bishop et al (1986) and Lundgren et al (1989), when the flexible theory 
provides an insight into unexpected happenings which would otherwise not be accessible. 


Slamming 


51 So far only loads occurring at wave encounter frequency have been discussed. However, 
waves can also cause loads at much higher frequencies due to impacts between the ship's hull and 
the water surface, known as slamming. Clearly, this type of response will usually occur when 
some part of the bottom of the ship comes out of the water and then re-enters, but similar 
responses can be seen resulting from very rapid changes in added mass and damping due '° 
sudden variations in hull velocity and acceleration. Slamming in theory can happen at any point 
along the length, but the most critical area will be at the keel in the forward part of the hull. 
Important slams can also occur at the stern in following seas, and at flare in the bow. The theory 


| Gus. el 
and means of load prediction are much the same in each case, so only the most severe case of ke 
Slamming forward will be covered here. 


. . tin 
52 | The immediate effect of a slam is a short duration high pressure pulse on the hit , ae 
Typically, at the keel, this pressure may rise to as much as 2 MPa in 0.05-0.1 sec, gue will 
further 0.2 sec. Where relative velocities between water and structure are lower the pres?" "oy 


; ‘ , 3e8 a 
be less but the rise time will be greater, up to 1.5 sec. This pressure pulse causes enone 
three structural levels, primary, secondary and tertiary. The primary response consists oF ™ icant 


vibration of the hull as a beam in all modes of which the lowest three or four will be sign 
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value. Originally the threshold velocity, Ss., Was stated by Ochi to be equal to 0.093 (Lg), but this 
does not hold for later experimental evidence and there is some doubt as to the real value for 
different ship types. 


55 If the two events noted above are denoted U and V and the draught at the forefoot is Re 
then the combined event (U 7 V), assuming they are independent has a probability 


(3.36) P(U A V) = P(U) .P(V) 
The probability of a slam is therefore 


(3.37) P(slam) = P(U > T;) .P(V > Sx) 


In the same way as response amplitude operators can be deduced for motions and loads, RAOs 
can also be found from strip theory or model tests for the two variables T, and s., and from these 
may be deduced the relevant response spectra. It follows therefore that 


ee) P(slam) = exp {ie 4 mal 


2M5, 2M, 


| 
5 
4 
4 


noting that m,, is the second moment of the response spectrum because S.. is a velocity. From this 
can be deduced the expected number of slams in unit time (N.) which is 


(3.39) N, = (1/{2n}) Vim,,/m3,) .P(slam) 


56 Current empirical methods only allow for an estimate of peak pressure at some point on 
the keel on normal hull forms. Ochi and Motter (1973), from a large number of wedge drop tests, 
postulate that the pressure p = ks, where k is a hull form parameter constant for a given section. 
As the probability distribution of s is known then the probability density function of p can be 
derived, hence the short term probability of p exceeding any limiting magnitude. 


57 Ochi’s conclusions from his drop tests are: 


a. k is independent of ship speed up to 15 knots for a 150 m ship. 
b. k is independent of sea state. 


G7 KI5%4 function of hull form, particularly Shape up to 10% of draught from base. 


d. A ‘V’ form is better than a ‘U’ form for pressure magnitude. 
e, Peak pressure moves aft with increasing speed. 
It should be noted, however, : 
ate ploted to sttaight liner on og ese the relationship in paragraph 56 and all his resus 
as ines on a log scale, and that all the results from ship model tests are lower 
48 
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than those from 2-D wedges at the same velocity. Further prediction of values of k by Ochi 
tewison and Belkin may be found in the report of ISSC Committee 8 ( 1973) | 


sg The time history of the pressure pulse may be ex 
both of which may be input to strip theory, although a flexible beam method is necessary due to 
higher frequency responses being significant. The first approach is known as impact slamming 
and is due to the short duration, high pressure pulse discussed by Ochi and Motter (1973) and 
refined by Stavovy and Chuang (1976). The second approach is known as momentum transfer 
slamming and covers a longer period (~1 sec) lower pressure load as water is moved outwards due 
to the downward motion of the ship; this phenomenon is discussed in detail by Leibowitz (1963). 
Both methods have been incorporated into computer simulations and are compared by Belik and 
Price (1982) and by Clarke (1986). The computer simulations have been used by Bishop, Clarke 
and Price (1984) to compare with results of comparative seakeeping trials between two frigates 
(see Clarke (1982)) and have been shown to give reasonable agreement, but the complexity of the 
input data required is such that the method is not yet suitable for use during the design process. 
Clarke (1986) has, however, shown that the effects of slamming on vertical bending moments 
extend significantly forward of amidships, and it is recommended therefore that the design 
bending moment lengthwise distribution, as derived either by static balance or by strip theory, is 
extended forward at a constant value from the maximum at amidships to a point at 0.35L from the 
bow, and then reduced linearly to zero at (or near) the fore perpendicular (see also Chapter 6). 


plored in two different ways, either or 


59 With regard to local response, Jones (1973) has shown that in general for monohull type 
scantlings it is permissible to treat the problem as quasi-static, as the natural frequency of 
clamped plates is much greater than the effective pressure pulse frequency. Consequently, for flat 
plates any standard analysis method may be employed, although it is probably best to assume that 
the plate edges are pinned but free to pull in, as the pressure at any given time will be localised and 
likely to be acting on only one plate panel. For the response of a grillage the same assumptions on 
frequency of response will apply, that is a quasi-static solution should be adequate and, because of 
the localisation of the load, it will usually be sufficient to apply the pressure to each stiffener and 
associated plating independently of other structure. A simplified method for estimating slamming 


Pressures is given in Chapter 10, together with methods of designing the bottom structure to resist 
the pressure 


60 
other p 
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The last aspect of hull response to slamming is whipping, which may also be caused by 
henomena such as underwater explosions (see also Chapter 4 and Chapter 15). As stated 
Whipping is the generalised response of the hull, vibrating as a beam, to any impulsive 
8 which in the current context is slamming. Dynamic analysis of non-uniform beams is now 
a fairly Standard process, for example Myklestad (1944) for aircraft wings, developed by Todd 
(1961) for Ship applications. It is therefore possible without undue difficulty to estimate 

“quencies and mode shapes for the dry hull, assuming that sufficient information is available on 
a ‘, ertia and stiffness distributions. However, extending that analysis to provide reliable 

Sults from impulsive loading while at the same time making proper allowance for fluid added 
te ; as damping is still, for practical purposes, unreliable although the theory exists, for 
| “ample. \shop et al (1978). (This is the same theory as the computer simulation mentioned 
Seq eae ‘ragraph 58.) It is still necessary to correlate such analyses with measurements Sie 
Sure that unexpected deviations in the results have not occurred in the numerica 
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Green Sea Loading 


61 Another effect of slamming, and indeed of proceeding at speed into even moderate seas. js 
to throw green water onto the deck of a ship. The worst affected parts of the ship will clearly be 
the forecastle and bridge front where quite high pressures and consequential damage may be 
experienced. However, green seas can affect any exposed part of the structure and pressures must 
be allowed for in design. It should be noted also that these pressures may occur at the same time ay 
those due to wind (Chapter 4) as well as wave bending loads and therefore combined load effects 
must be taken into account. 


62 The numerical treatment of pressures due to green seas is non-existent (although 
frequencies may be computed with some degree of accuracy on the same lines as slamming 
frequency) and so design pressures are derived from a knowledge of structures that have not 
failed, and the pressures that would have been needed to fail them. A few measurements have 
been taken at sea but these are mainly using small pressure transducers which measure only very 
localised peak pressures and are meaningless for design of structure. Consequently there is 4 wide 
variation in design pressures quoted by design authorities and classification societies and the 
requirements of the latter, for particular parts of ships, may differ by an order of magnitude. 
Traditionally for RN warships a figure of 50 kPa (100 Ibf/ft2) has been used for lack of better 
information, and this figure finds some support in the work of Yagle (1968) who studied the 
pressure on decks, side shell and deckhouses on a destroyer model in irregular seas. At 00 time 
was a pressure greater than 50 kPa measured. In general, this figure may be too severe, but there 
some evidence from recent damage on British frigates that it may be too low for the bridge fron! 
and for equipment on the forecastle. The following design figures are therefore recommended for 
frigate and destroyer type vessels; for smaller and larger ships recourse should be had a 
Classification Society rules for vessels of about the correct size and speed. All pressures should 
used as if they are applied statically over an appropriate area of the structure, for example pa 
_ any 25% of the bridge front. On decks, while a significant area of plating may be subject sitet 
pressure at one time, only a few stiffeners will be loaded by any single green se. It is pre 
adequate, therefore, to halve the load when designing stiffeners because the plating 
average out the load. | vs : 


| 4_Bridge front and forecastle mounted equipment - 100 kPa (but see 8) 
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Upper deck forward of the bridge front - 25 kPa at the bridge front, rising linearly to 


1.0(VL - T,) kPa (if it is greater than 25 kPa) at the extreme fore end where L (LWL) 
and T (draught at forefoot) are in metres. 


c, Side shell between waterline and upper deck - 50 kPa. 


d. Exposed areas of the upper deck and 01 deck - 25 kPa. 


e. Exposed decks from 02 upwards - 15 kPa. 


Exposed vertical surfaces other than the bridge front - to vary linearly from the 


pressure required for the contiguous side shell or deck at the bottom to that required 
for the deck over. 


ge. High value equipment, for example weapons in exposed locations such as the 


forecastle, may need to be designed to withstand up to 200 kPa to minimise the risk of 
loss or unacceptable degradation. 


Horizontal Bending and Torsion 


64 Both horizontal bending and torsion are caused by wave action when seas off the bow are 
being met. Both can be given the same type of analytical treatment as vertical bending except that 
the anti-symmetric degrees of freedom will be needed. An example of use of the theory for these 
responses, albeit for higher frequency distortional modes, is given by Bishop et al (1986). 
However, for conventional monohull vessels, stresses due to both these responses are much lower 
than for vertical bending and can usually be ignored, especially as their maxima are likely to 
occur when the vertical response 1s well off maximum, although it has been shown on a trial that 
some vertical and lateral responses 1n moderate seas can be in phase with a ratio 2.5:1. An 
exception may need to be made for ships with very large deck openings (for example container 
ships and possibly warships, ‘f the machinery downtakes and uptakes are particularly large) and 
for ships of very shallow section when torsional stresses may be significant. 


i little available information on which to base estimates of wave induced 
feos e abe hae recent is that suggested by Meek et al (1972) who propose for the OCL 
container ships a figure of 0.0055LBD? tonne.metres maximurn at amidships with a (I -cosine) 
distribution lengthwise. however, there is some evidence from strip theory that rolling has a 
Significant effect on torsional loads and that the maximum could be reduced slightly but should be 
extended some distance either side of amidships producing a trapezoidal form of distribution. A 
good compromise would appear to be a maximum torque of > x 10-5 LBD? MNm over the centre 
0.6L tapering linearly to zero at the ends, but further discussion and other proposed formulae are 
presented by Vedeler (1924) and De Wilde (1967) together with some experimental verification 
by Lloyd et al (1980). 
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Figure 3.10 Hydrostatic Heads Recommended for the Design of the Bottom and Sides 


Conclusions 


7) Anextensive literature, both historical and contemporary, exists in support of quasi-static 
estimates of design wave bending moments. The majority of these, however, relate to some 
arbitrary severe, but not extreme, condition which fits neither a serviceability nor an ultimate limit 
state. The ultimate limit state which is now recommended for design and which will compare with 
ultimate strength estimates to follow (Chapter 7) is covered by an update of the quasi-static 
estimate which allows implicitly for slamming and other potential inaccuracies by means of 
correlation with measurements from sea and with corresponding partial safety factors to pracesie 
a design bending moment with a known probability of exceedance in a given number of wave 

This approach allows for non-linearities in both ship and wave form but can only 
romances P ‘c effects of motion implicitly through the correlation with measurements. 
Consequently, it cannot take account of novel or unusual hull forms and, while different 


operational scenarios can be allowed for in the mathematical model, any extrapolation from the 


‘existing database must be treated with great caution. 


deal of effort has been put into dynamic analysis of ship 
for motion and for loads, and allowing for the — 
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